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I. INTRODUCTION 



As already stated in many publications before, grav- 
itational waves from inspiralling compact binaries are 
the most promising sources for ground based planned 
and already operating gravitational wave (GW) detec- 
tors. To guarantee successful search for GWs, one needs 
to obtain promising search templates incorporating all 
important physical effects that have an influence on the 
form of the signal. Ground-based detector networks 
like LIGO (USA), VIRGO (France/Italy) and GEO 600 
( Germany/UK) have the sensitivity to be able to see the 
last seconds or minutes of the binary's inspiral, where 
the corrections, coming from general relativity, to the 
Newtonian orbital motion get important, depending on 
their masses. In order to detect GWs from inspiralling 
compact binaries without spin in quasi-circular orbits, a 
large library of ready-to-use inspiral templates has been 
put up pp. Eccentric inspiral models without spin have 
also been developed [2J [3J and are well understood. Re- 
cently, Yunes and collaborators obtained a formalism for 
frequency domain GW filters for eccentric binaries [4|. 
All of them heavily employ the post Newtonian (PN) ap- 
proximation to general relativity. 

It has been shown by several authors, that for a suc- 
cessful detection, effects of spin have to be included and 
foundations for the detection of spins have been laid down 
[5J EJ 13 IS] • During the inspiral phase, before reaching 
the last stable orbit, those effects are long-term modu- 
lations of the GW signal in a comparison with the time 
scale of only one orbit. They can lead to substantially dif- 
ferent shapes of the signal compared to those ones show- 
ing up if the spins are neglected. The foundations for 
the motion of spins in curved spacetimes are given in 
[Ttj] . In harmonic coordinates, the spin dependent EOM 
were derived up to next-to-leading order in the spin-orbit 
coupling by Faye et al. pj] and Blanchet et al. [12] . 
where velocities have been used to characterize the or- 
bits. In Arnowitt-Deser-Misner coordinates [13], higher 



order Hamiltonians dictating the equations of motion for 
orbits and spin (from this point on referred to as EOM) 
have recently been derived by Damour et al. |14j and 
Steinhoff et al. [TH1 EB] . The spin-independent part of 
the binary Hamiltonian is known to 3PN order [IT] . 
The solution to "simple precession" of the leading or- 
der spin-orbit interaction, which was the case for single 
spin or equal mass, was discussed in |9 and later in |18j . 
where the GW polarizations h x and h + were derived as 
a PN accurate analytic solution for eccentric orbits. The 
latter has heavily inspired this work, which will give an 
approach to the more general case of unequal masses and 
arbitrary two-spin configurations. 

The paper will be organized as follows. Section [TT] will 
present the involved Hamiltonians and the associated 
EOM for the binary in the center-of-mass frame. In 
section |III| the geometry and the coordinates relating 
the generic reference frame with the orientation of the 
spins and the angular momentum vector are provided and 
characterized by rotation matrices. The time derivatives 
of these rotation matrices will be compared by Poisson 
brackets in section ITVl and first order time derivatives of 
the associated rotation angles will be obtained. A first- 
order perturbative solution to the EOM for the spins 
is worked out in section [V] The orbital motion will be 
computed, for quasi-circular orbits (circular orbits in the 
precessing orbital plane), in section VI As an applica- 
tion, the resulting GW polarizations, h x and h + in the 
quadrupolar restriction, are given in section [VII| 



II. THE CONSERVATIVE EQUATIONS OF 
MOTION FOR THE SPINS 



In this section the dynamics of spinning compact bina- 
ries is investigated where the spin contributions are re- 
stricted to the leading order gravitational coupling. The 
Hamiltonian associated therewith reads 



1PN 



H 



2PN 



T^so , 



(2.1) 
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with 7Yn, Wipn and H2PN respectively are the Newto- 
nian, first and second PN order contributions to the con- 
servative point particle dynamics (e.g., [19 j and refer- 



2 



ences therein) and Tiso is the leading order spin-orbit 
Hamiltonian [20] . 

In the following computations, use will be made of the 
following scalings to convert the quantities in calligraphic 
letters to dimcnsionlcss ones on the rhs: 
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(2.4) 

(2.5) 
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where m a is the mass of the a th object (a = 1,2), m is 
the total mass, m = mi + m 2 , A* is the reduced mass de- 
fined as mi m 2 /m and the symmetric mass ratio is given 
by 77 := mi m 2 /m 2 . The variables p and r are the scaled 
linear canonical momentum and position vectors, respec- 
tively, and commute with the spins S a - Explicitly, the 
contributions to the scaled version of Eq. (2.1 ) read 



H(r,p, S U S 2 ) 



with 



H N (r,p) + H 1PN (r,p) 
+ H SO {r,p, Si,S 2 ), 



H 2PN (r,p) 
(2.6) 



(1 - 5?7 + brj 2 ) (p 2 ) + - (5 - 20r? - 3t? 2 ) (p 2 ) - 2n 2 (n ■ pfp 2 - 3?7 2 (n • p) L 
8 L 



(2.7a) 

(2.7b) 
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r 

(2.7c) 
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where r = \r\ and S c f[ is the so-called effective spin, 



S e s = SiSi + 5 2 S 2 
3 



61 = 



1 - y/1 - 477 



(2.8a) 
(2.8b) 

(2.8c) 



Considering only the spin-independent part of the Hamil- 
tonian, the orbital angular momentum vector is a con- 
served quantity. The motion of the reduced mass /i will, 
without SO interactions, take place in a plane that is 
perpendicular to L and that is invariant in time. Adding 
the spin-orbit term will, in general, lead to a precession 
of the orbital angular momentum. The EOM for L, de- 
fined by L := r x p and the individual spins Si & S2 
can be deduced from the equations 



— = {L, H so } — -*-^Sett x L , 
at c z r 6 

= {Si,H so } = ^3-^ x ^1 j 



dt 
dS 2 



{S 2 ,H^o} — ~b~^^ x ^2 



(2.9a) 
(2.9b) 
(2.9c) 



Equation (2.9a I describes the precession of L w.r.t. 
the total angular momentum vector </, defined as 
J = L + Si + S 2 . The key idea in the next sections is 
to compute time dependent rotation matrices for L, Si 



and £2 for a number of rotation axes and angles that are 
to be introduced in the next section. Let us state that 
the magnitudes L, Si and S 2 of the vectors L, Si and 
S 2 are conserved, 



dL 2 _ 
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~dT ~ 


dt 


dS 2 _ 
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(Si -Si) 
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2So 



L-(S eS xL) = 0, (2.10a) 



Si-(Lx Si) = 0, (2.10b) 



S 2 -(Lx S 2 ) = 0. (2.10c) 



Equations (2.9) show that L 



(Si + S 2 ) and, thus, 
(2.11) 



the total angular momentum vector J satisfies 

dJ d\J\ 

- = 0, gl vmg — = 0. 

The magnitudes of S and S e g behave as follows 
dS 2 3 % /T^4^ 



dt 

dS^s 
dt 



. , L-(S 1 xS 2 ), 

c z r a 

3yi~^4^(12 + ?7)7/ 



4c 2 r 3 



(2.12a) 

L-(SixS 2 ). (2.12b) 



Notice the conservation of S 2 S in both the test-mass 
(7/ = 0) and equal-mass (77 = 1/4) cases. Using above 
equations, we will be able to compute the evolution equa- 
tions for the rotation angles. The associated geometry is 
introduced next. 



III. GEOMETRY OF THE BINARY 

As done in [18], it is very useful to use a fixed or- 
thonormal frame (ex,ey,ez) and to set &z along the 
fixed vector J. The invariable plane perpendicular to 
J will then be spanned by the vectors (ex,ey). The 
motion of the reduced mass will take place in the or- 
bital plane perpendicular to the unit vector k := L/L. 
For a clear understanding of the following, please take a 
look at Fig. [I] First, the vector k is inclined to ez by 
the (time- dependent) angle O, which was also the con- 
stant precession cone of L around J for the single-spin 
and equal-mass case of [13]. As before, the orbital plane, 
itself spanned by the vectors where j = k x i, in- 

tersects the invariable plane at the line of nodes i, with 
the longitude T measured in the invariable plane from 
ex- 

The geometry of the binary will be completed by the 
spin related coordinate system (i s ,j s ,k s ). This frame 
is constructed from the system (i,j,k) to be rotated 
around the axis i to point from the top of L to the top 
of J with the new direction k s . In other words, this spin 
coordinate system is chosen in such a way that the total 
spin, Si + S2, has only a k s component and i s = i holds. 
If is known, the spins are left with an additional free- 
dom to rotate around k s by an angle (j) s (the index "s" 
is a hint for positions in the spin system). This angle is 
measured from i a to the projection of Si to the (i s ,ja) 
plane, similar to T's function in the reference frame. 

There exist simple geometrical relations that will re- 
duce the freedom to choose rotation angles arbitrarily, as 



will be shown in the next subsection. 



A. Geometrical issues 

As mentioned already, in this geometry the spins and 
angular momenta - being fixed in their magnitudes - 
only have three degrees of freedom: the angles @, T and 
(f> s . Once © is determined, also a^ s (the angle between L 
and S) is fixed and so is magnitude S of S = Si + S2 by 
triangular relations. Calling an the angle between the 
spins Si and S2, the following equations list the rotation 
angles and magnitudes as functions of , where also use 
is made of the sin relations, 

5(0)= VJ 2 -2JLCOS0 + L 2 , (3.1a) 
a 12 (0) = cos- 1 ( g(e) ^ 2 ~^ ) . (3-lb) 

a ks (0)= .-sin-^^M) , (3.1c) 

-vr>s, • -1 /S , 2 sina 12 (0)\ 

s(Q) = sm { m j . (3.1d) 

These relations will be used extensively to simplify the 
angles evolution equations. How they are incorporated 
and applied will be shown next. 



B. Coordinate bases and associated transformation 
matrixes 

This section introduces the coordinate transformations 
from the reference system to the orbital triad and the spin 
system. To construct the EOM for the 3 physical angles 
0, T and <J) S , the idea is to compare the evolution of these 
rotation angles - as arguments for rotation matrices - 
with the Poisson brackets, Eqs. (2.9a) - (2.9c I. Let us 



begin with the explicit computation of the transformed 
coordinate bases. 

1. The orbital triad (i,j,k) can be, not surprisingly, 
constructed by only 2 rotations from the reference 
system. In terms of rotation matrices, we have 




1 
cos 
-sin© 




(3.2) 





The spin system is constructed, simply by another 
rotation of a^ s around the vector i, from the orbital 
triad, 



(3.3) 



such that i s = i holds. Important note: the angle 
ctks has negative sign relative to 0. That's because 
k s has to be moved "backwards" to point to J! 

Having transformed the unit vectors with these ma- 
trices, the coordinates transform by their transposed in- 
verses, which are - in case of rotations - the matrices 
themselves. 

Now, we have everything under control to construct 
the set of all the physical vectors. I will list all of them 
below. First, let me define some shorthands for rotation 
matrices: 



'10 
[0] = (0 cos© sin© 
,0 — sin cos © , 



(3.4) 



FIG. 1: Binary geometry completed by a rotating spin coordinate system. The usual reference frame is (e^,ey,ez) having 
chosen ez to be aligned with J. The vectors L,Si,S2 describe the orbital angular momentum and the individual spins, 
respectively. The angle Q denotes the inclination angle of L w.r.t. J, which is - of course - to be taken as a time dependent 
quantity. The orbital plane, being perpendicular to L by construction, is spanned by the orthonormal vectors j and i, where 
the latter one intersects the invariable plane at the angle T measured from ex- The spin-coordinate system is constructed 
out of the orbital dreibein (i,j,k) by a rotation of a ks around i, such that the vector pointing from I to J is the total spin 
Si + S2- The angle ai2 is measured between Si and S2. The spin Si, projected into the (j s ,i s = i) plane is rotated by an 
angle <f> a from i, and Si itself is moving on the circle (with variable radius) embedded in the figure. 



cos T sin T 0\ 
[T] = ( -sinT cosT , (3.5) 
1/ 

fl Q \ 

[a ks ] = cosa ks -sina ks . (3.6) 
y sin a ks cos a ks J 

The orbital angular momentum L in the reference system 
(indices labeled inv) arises from two rotations from the 
orbital triad (ot) where it has only one component: 

L = {[6(t)] [Tit)}}- 1 (0,0, L), (3.7) 
or, in components, 



(£)r= [6(f)] [T(i)] (L) 



(3.8) 



The spins, in the spin system (s), where the fc s is aligned 
with S = Si + S2, have the following form, 



Si = Si (cos 4> s sin s i s + sin <p s sin s j s + cos sk s ) , 

(3.9a) 

S 2 = Sks-Sx, (3.9b) 
S= Sk s . (3.9c) 



IV. THE TIME DERIVATIVES OF T, 9 AND 4> s 

To obtain an EOM for the angle O, one possi bility i s 
to use the time derivative of \S\ 2 = (S ■ S), Eq. ( 2.12a l, 
to apply this, for example, in the spin system and to 



compare the result with the time derivative of Eq. (3. la I 
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with = Q(t). The result is The same result will be obtained by computing the time 

derivative of the orbital angular momentum L in the in- 
£ = -_---? sin a ks cos ^ s sins esc (4.1) variable system. Therefore, take Eq. computers 
2 J time derivative and finally compare the result with (2.9a I . 

. , Because the angular velocities appear in relatively simple 

relations, it is easy to extract them from the ex and ey 
_ An entry. The results are 

^s = -^A (4-2) 



T = — C^cscG [5i(<$2 — <^i) cos aks sin <p s sin s + sin aks(Si(<5i — 82) cos s + S62)] , 
= Cl Si(5i — 82) cos(f) s sin s , 



(4.3) 
(4.4) 



with 
of a ks 

metrical relations, Eqs. (|3.1 



Cl '■— (c 2 r A )~ . The functional dependencies 
. s and S on O are implicated. Inserting the geo- 
it turns out that Eqs. (4.1) 



and (4.4 1 are equivalent. Also, the allegedly worrying 



asymmetric appearance of the quantity Si can be stu- 
diously avoided by replacing s by its function of 0. 1 

Also note that, if the rel ation s r\ = 1/4 or Sj — (i — 
1 or 2) are inserted in Eq. (4.3 1, one recovers Eq. (4.32) 

of mo. 

Now, let us turn to the last quantity to be determined, 
the angle <f> 8 . The geometry offers various possibilities to 
calculate the time derivative of this angle. The easy way 
is to compute Si in the invariable system. In compo- 
nents, we have 



(Si)\ 



[aUt)} [©(*)] TO] (Si) 



(4.5) 



The time derivative of (4.5 I might be compared with Eq. 
(2.9b). The result will be given in terms of the angles 



already determined: since we already know and T on 
the one hand and aks as a function of on the other, we 
have the expression under full control. 

The other way is to take the Leibniz product rule for 
S, namely d t S = d t S l e, + S l d t &i with e t = (i a , j a) k s ). 
We we already know that i s = i, k a \\S and j s = i s x k 3 , 
whose time derivatives are already known. Both consid- 
erations result in 

4> s = C e (0 - a ks ) + C~ s a +n , (4.6a) 

Cq = tan <f) s cot(aks — 0) + sec(/>cot s , (4.6b) 

Cg = — sec 4> a cot (aks — 0) — tan <j> cots, (4.6c) 

n = - Ci L sin csc(a ks - 0) ■ (4.6d) 

For the case of equal masses (81 = 82 = 8 = 7/8), one 
obtains for 4> s , T and a very simple system of EOM, 

= 0, (4.7a) 



1 The angular velocities, | |4.3| l and | |4.4 1, are in complete agreement 
with Eqs. (5.11a) and (glib) of [21 . 



7J 

8c 2 r 3 ' 
LS sin 



esc < sin 



J sin© 







which can be integrated immediately, giving 



e(t) = e , 

T(t) = Q r t 

<j> a (t) = n$t- 

with the angular velocities 



To, 

0s O 1 



L(5sin0 o i • -i 
llrh = s— 5 — csc < sin 



J sin 0n 



S(Q ) 



0o 



(4.7b) 
(4.7c) 



(4.8a) 
(4.8b) 
(4.8c) 



(4.9a) 
(4.9b) 



Summarizing the EOM for the coordinate transformation 



angles, Eqs. (4.3), (4.4) and (4.6 1, this system of EOM 



can be written in a compact manner. Calling the vector 
of constants, C = {E, Si, 52, L, k ■ S e s} - where E and L 
are related in the case of quasi-circular orbits - and the 
vector of dynamic variables, associated with spins and 
angular momentum, X = {0, T, <fi s }, we may write 



dt 



X = Y C (X) . 



(4.10) 



A perturbative solution will be given in the next section. 



V. FIRST ORDER PERTURBATIVE SOLUTION 
TO THE EOM FOR THE NON-EQUAL MASS 
CASE 

The EOM for (0, T, (j> s ) can also be solved by a simple 
reduction scheme. We assume that the deviation from 
the equal-mass case is small compared to unity, 



Si 



Si 



< 1 



(5.1) 
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Then, having the equal-mass case under full analytic con- 
trol, we can construct a perturbative solution to the non- 
equal mass case. The proceeding is as follows: Imagine a 
system of EOM for a number N of dependent variables 
X: 



X = Y(X) . 



(5.2) 



The time domain solution to this system is denoted by 
the superscript "0" , viz 



X(t) = X (0) (i) 



(5.3) 



Let us assume that the EOM, Eq. ( |5.2| , are perturbed by 
some terms of the order e (e is a dimensionless ordering 
parameter) , 



X = Y(X) + eP(X) 



(5.4) 



The solution at the first order in e can be obtained by 
adding a small perturbed quantity to be determined to 
the solution of the homogeneous equation, 



xP(t)=Xr(t)+eS i (t). 



r(0) 



(5.5) 



Inserting this into Eq. (5.4 1, one obtains 



xP 



x (0) 



eSi 



Yiixf + eSA + ePiiXf'+eSj) 



(o) 



€ ^dx~ 

3=1 J 



+eP l {xf ) ) + 0{e 2 ) (5.6) 
Comparing the coefficients of the two orders of e gives 
0: Xf = Y{xf ] ), (5.7) 



N 



3 = 1 



dXi 



(5.8) 



The first equation is solved via definition, and what re- 
mains is the second, having inserted the unperturbed so- 
lution in the perturbing function P. For our purposes, 
N = 3 with X = {T, 6, <j> s } is a small number of EOMs, 
but complicated functional dependencies are included. 
The matrix appearing in Eq. ( 5.8 ) does not mean a prob- 



lem to us, because fortunately, the only dependency of 
the sources is on O. 

For our computation, we need to divide the EOM into 
a non-perturbative and a perturbative part. In the fol- 
lowing, we use the definitions 



Xi 



X2 



s 1 + s 2 

2 

Si - s 2 



(5.9) 
(5.10) 



Rewriting the EOM for the angles in terms of xi and %2, 
labeling all \i contributions with the order parameter e 
as well as inserting the non-perturbative solution, Eqs. 
(4.8) to these terms, one obtains 



9(i) 

j(0) 



e C L Si2\2 cos(t fl,/, + 4> ) sin s(6 ) , 

ClS(Q)xi sin ak s (©) esc O +e [ClX2 csc 6q (25i cos aks (©o) sin s (0q) sin (tfl^ 



=Cl Jxi— const. 

sin a ks (9 ) (5(0 O ) - 25i cos s(0 o )))] , 

ds 



CiLsin0cscaks(0) + e Cg(Q,<f> s ) 
(Xi + e X2) sin csc a ks (0) + e 



dO 

C S ~(0,0 S ) 



Ce(0,^ s ) 
ds 



_ da ks 
dQ 

C e (0,0 s ) 



_ da^s 

oe 



e. 



e=e 



(5.11a) 

(5.11b) 
(5.11c) 

0. 

(5.11d) 



The parameter e simply counts the order of the perturbative contribution and is later set to one. The first term for 
T is constant and thus does not have to be expanded in powers of e, but the associated first term for <p s does, such 
that the perturbative solution for has to be included. Taylor expanding this term, removing all contributions to 
the unperturbed problem, what remains is a system of EOM for Sq, St, S4> that can be simply integrated, because 
as soon as S@(t) is known, all the other contributions are straightforwardly evaluated. Requiring that the perturbing 
solutions vanish at t — 0, the solutions are simply given by 



Se(t) 



SqcH . 



S r (t)= / Sydt 



(5.12a) 
(5.12b) 
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and explicitly read 



Se(t) 
Sr(t) 



S<p(t) = / S^dt . 
Jo 



— - — sin ai2 (0 ) ( sm <Psq - sin («2 + </> s0 )) , 



ClX2 csc6 



251 cos aks(o) si n S(o) ( cos ^so — cos + 

- if^ sin a ks ( ) (5(0) - 25i cos S( ) 
S*(t) = C stat * + C (t) + Ci(t) + C e (t) + C a (t) 



$30)) 



with the shorthands 



CoM = 



C g (t) 



Xi 



C'lJSiS 2 X2 sin6 sinai2(o) (ifi^sin.0 ao + cos (tf2^ + <fi s0 ) - cos<fo) (-2c 2 Jr 3 fT - 2xiS 



xiSf 0) n^sf 0) -j^m 2 e 

-2C L Jx2 cotai2 (0 ) - S , iS , 2 sinQ;i2( 



2(0) 



XiSf 0) fl<j> J Sf Q) - Sf sin 2 ai2 (0 ) 



c (t) 



(c 2 r 3 tr20 sin S(o)^* + £xi sm ©o cos S(o) (cos c/) s0 — cos (tQ^ + </> s q))) > 



2S\tx2 cos S( ) 251X2^* esc 6 sin 5( ) (cos SO - cos (ifi^ + </> s0 )) 



Lxi£l 4 



C a {t) = C e (t) 



J ( Sis cos 6 - JX sin 2 6 



o2 
*(0) 



Sf 0) - J 2 sin 2 e 



the initial values of the functions (3.1al - ( 3. Id I 



S( ) 


= 5(0 O ) , 


aks(o) 


= aks(B ) , 


"12(0) 


= ai 2 (6o) , 


S(o) 


= ~s(&o) , 



and the definition 



VI. THE ORBITAL MOTION 

The motion of the spins is only half of the physical con- 
tent of the spin-orbit dynamics. Once we fully have the 



L 2 Xl sin 2 e 



c 4 r 6 ft 2 



r 



(5.12c) 
(5.13a) 

(5.13b) 
(5.13c) 

(5.14a) 
(5.14b) 



(5.14c) 



(5.14d) 



(5.15a) 
(5.15b) 
(5.15c) 
(5.15d) 
(5.15e) 



(5.15f) 



motion of all the spin-related angles under control, we 
might turn to the orbital dynamics, i.e. the motion of 
the reduced mass in the orbital plane. It will turn out 
that employing coordinate transformations will be very 
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helpful here, too. 

The aim is to solve the orbital EOM to the full Hamil- 
tonian, 



H = H* 



Hi 



PN 



H- 



2PN 



so 



(6.1) 



At this point, we can do a useful simplification. As long 
as we incorporate only leading order spin dynamics, only 
Newtonian point particle and spin dependent contribu- 
tions will mix at the end, higher order PN terms cou- 
pling with the spins will be neglected consequently. For 
the computation of the spin dependent part of the orbital 
phase, therefore, we only have to take -ffN,so = H^+H^o 
and add the 1PN and 2PN (spinless) terms for the point 
particle afterwards. 



H — i?N,SO + -ffiPN + #2PN 
V = <y5N,SO + V?1PN + ¥>2PN • 



(6.2) 
(6.3) 



The Newtonian and spin orbit part of eq. (6.1 1 reads 

,2 1 



H 



N,SO 



El 

2 



1 



(rxp) ■ S 



off • 



(6.4) 



and can be handled with the method described in [15] , 
The aim there was to introduce advantageous spher- 
ical coordinates, (r,9,<fi), with their associated ONS 
(n, eg, e^) with e z • n — cos 9, n ■ e x = cos <j> sin 9, as 
can be seen in Fig. First, we define the normalized 
relative separation vector according to 



n = sm 9 cos <p ex + sin 9 sin <p ey + cos 9ez 



(6.5) 



The time derivative of r, the linear momentum p, its de- 
composition in radial components and the corresponding 
orthogonal ones can be written as 



r = rn , 

r = rn + r9 eg + r sin 



p = p r n + p g e s -t 

2 i „2 



1 2 
P =Pr 



I? 



(6.6a) 

■4 , (6.6b) 
fie<f, , (6.6c) 
(n -p) 2 + (n x pf 

(6.6d) 



Inserting p into Eq. (6.4 1, computing p^ — p ■ e^ and 



using the orthogonality relation of the used triad, one 
obtains 



P<f, 



2E- 



I? 



2(L ■ S cff ) 



r sin 9 

L 2 2 

„2 v 4> 



L 2 - 



sm 



2, 



(6.7a) 
(6.7b) 
(6.7c) 



In [T5] , it was possible to reduce these equations by some 
algebraic relations and the fact that the angle was 



J = Je z 




L = Lk 



FIG. 2: The geometry of the binary, having added the ob- 
server related frame (p,q,N) (in dashed and dotted lines) 
with N as the line-of-sight vector, after removing the angles 
in the spin frame. The line-of-sight vector is chosen to lie in 
the ey — ez-plane, and measures an angle io (associated with 
the rotation around ex) from ez, such that p — ex, and this 
is the point where the orbital plane meets the plane of the 
sky. Because of this rotation, the angle io is also found be- 
tween the vector q, itself positioned in (ey , ez), too, and ey . 
The grey area in the graphics completely lies in the orbital 
plane, spanned by (i,j) and <p measures the angle between 
the separation vector r and i. The polarization vectors p and 
q span the plane of the sky. The inclination of this plane with 
respect to the orbital plane is the orbital inclination i. The 
inclination of the orbital plane with respect to the invariable 
plane is denoted by O. Please note that L does not lie on the 
unit sphere, only k does! 



constant in time - here, it is more complicated. It is still 
allowed to express L z , the projection of L onto ez, in 
Eq. (6.7b I and (6.7c I over with the help of 



L cos 



r sin 9 

L2 



1 



COS 2 

sin 2 9 



(6.8a) 
(6.8b) 



Above equations are, for our purposes, the most simpli- 
fied versions of the p components and will enter in the 
dynamics of the angle ip in their current form. 

Our aim is now to connect the coordinate velocities, 
namely r, 4> and 9, to conserved quantities associated 
with the Hamiltonian of Eq. 



(|6_4 

velocity in spherical coordinates, 



ing formulae using Hamilton's EOM, r 



Computation of the 
5q. (6.5 1, gives folio w- 
dH NSO /dp, 



n x eg = e^ and n x e 



rsm( 



r = n ■ r = p r , 

r9 = eg ■ r = pg - 

>0 = • r = 



Setf 



eg ■ S. 



off 



(6.9a) 
(6.9b) 

(6.9c) 
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Of course, in the case of quasi-circular motion, f = — p r 
holds for all times. Remembering the geometry of Fig. 
[2] we recall that r is lying in the plane orthogonal to L, 
which itself is spanned by the vectors i and j. Calling ip 
(the orbital phase) the measure for the angular distance 
from i, we can write 



r = r cos ipi + r sin <p> j 



(6.10) 



The comparison of r, given by Eqs. (6.6a I and (6.5 , with 



the one in the new angular variables, Eq. (6.10) with 
Eqs. (3.2), implies the transformation 



where S$ is a shorthand for S c s ■ e^. The ambiguity of 
the sign in the first term can be removed if one takes the 
rotation sense of the reduced mass, or equivalently, the 
direction of L into account. Having (initially) the vector 
L in the northern hemisphere, one should choose "+" in 
above equation. This condition then holds anytime as 
long as Si + S 2 < VL 2 + J 2 . 

The quantity L/r 2 represents only the Newtonian 
point particle contribution. To express r and L in Eq. 
(6.131 in terms of E, one only needs Newtonian order, 



(M)-(T,¥>): 



cos 9 = sin if sin © 
sva.(4> — Y) sin 9 = sin ip cos © 
cos(</> — T) sin 9 = cos ip . 

(6.11) 

Time derivation of the first equation will give an expres- 
sion for 9, which can be simplified using the third one. 
The final expression is 



sin A 6 



cos 2 6 . 

1 -^tt:^ 



(6.12) 



with A = (f>—T. Setting O constant, one naturally recov- 
ers E q. (4 .28a) of [18]. Using this equatio n to eliminate 
9 in (6.9b I and after substition ±pg from (6.8b), one ob- 



tains a solution for (p and T 



S,i, 



sin A 



1 - 



cos 2 e c 2 r 3 

sin 2 6 



cos 2 e 
sin 2 6 



6, (6.13) 



r 
L 



(-2E)- 1 , 
{-2E)- 1 ' 2 . 



(6.14) 
(6.15) 



Summarizing the evolution for ip, one can separate it into 
a pure point particle (PP) and the spin orbit part (SO), 



(f = (fpp + (p so 



(6.16) 



The full 2PN expression for ippp can be extracted from 
Eqs. (5. 6c), (5.6d) and (5.6k) of [IB] without spin depen- 
dent terms, 



n=(-2E)^\l+ ^(-15 + ^) + ^ 

35 - 1077 



c 2 L 2 
setting e t = in 



Ac 2 



-2EL 2 



e 2 = 1 + 2EL 2 



-2E 



555 + 30?/+ llr? 2 



192 



V-2EL 2 



(5 - 277) 



877 + (17 - 7r])(-2EL 2 ) - 8xso cos . 



S 



^ j 8 + 4r; + 20?7 2 - (-2£;i 2 )(112 - 47?7 + I677 2 ) + 24y / -2EL 2 (5 - 2rj) 



(I7-II77) 



24 



(-2EL 2 ) v ,; V-2EL 2 

to eliminate L and using ippp = n (1 + k) |22j . giving 



(5 - 277) 



ipPP = {-2Ef' 2 { 1 + 5 (9 + 77) 



(-2E) 



<Pso 



-3(k ■ S eS )(-2E) 



3 , (-2S) 3 ^ 



891 201 11 2 
128 ~ "64~ 77 + 128 V 
sin AO 



(-2EY 



cos 2 e 



1 



cos 2 €) 



with 



(6.17a) 
(6.17b) 



(6.18) 

(6.19) 
(6.20) 
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cos(0 — T)[5i sin S sins(<5i — <5 2 ) cos(6 — a^ s ) — sin(6 — a^ s )(Si cos s(<$i — #2) + 5(0)52)] 

+ sin(</> — T)5i cos S sin s(<$2 — 5i) ■ (6-21) 



The first term in ipso, Eq. (6.20), comes from spin-orbit contributions to the value of L, as this is obtained from the 
energy expression (2.6), see section IV of [18]. The angle cf> can be computed with the help of Eq. (6.11) according to 

(6.22) 



4> = T + arccos ^cos cp/yX — sin ip sin 

Inserting the solutions 0(t),T(t) and 4> s {t) from sections IV and [V] to Eqs. (6.191 and (6.201, <p can be obtained by 
numerical integration, 



<p(t)= / <p dt + ip = (pp P t + / <p so (t) dt + <p . 
Jo Jo 



The radial separation at 2PN accuracy, after eliminating L, reads 



1 

F2£) 



(~2£) 
4c 2 



-7 + A(k-S cS )^(-2E) 



(~2^) 2 
16c 4 



-67 + 77(54 + 7?) 



(6.23) 



(6.24) 



VII. GRAVITATIONAL WAVEFORMS 

The gravitational wave polarization states, h + and h x , 
are usually given by 



1 



— (run; 4- n^r/A 1 



(7.1a) 
(7.1b) 



where and qi are the components of the vectors p and q 
orthogonal to the observer's direction, respectively, and 
hj^ is the transverse and traceless part of the radia- 



The leading order contribution, hJ T 



tion field , , 

where the subscript Q denotes quadrupolar approxima- 
tion, reads 123 



TT , 



4G/j 



GM 



(7.2) 



with Vkmij{N) as the usual transverse-traceless projec- 
tion orthogonal to the line-of-sight vector N, R' as the 
radial distance to the binary, the shorthands = Vi Vj 
and riij = Uj rij , using v = dr /dt as the velocity vec- 
tor and n = r/r as the normalized relative separation, 
respectively. 



Using Eq. (7.2 1, one may express both amplitudes of 
h x and h + as 



h+L = 



2G/j 

4 Ui 



i /',/', - qiQjj //, 



g{(P^-ta-) ! 



GM 



?• 

2G/x 



[(p • n) 2 - (q • n) 2 ] J , (7.3a) 



(PiQj +PjQi) ( v. 



GM 



c 4 ^ 



l -Up-v)(q-v) - ( ^(p-n)(q n) j . 

(7.3b) 



To compute the two gravitational wave polarizations, one 
requires an expression for the radial separation vector r 
and its first time derivative. It is efficient to give r ex- 
panded in the observer's triad (p, q, TV). In jTHj, this was 
done by expressing r in (ej,ey,ez) first, and secondly 
to compute this base from (p, q, N) as rotated around p 
with the (constant) angle Iq. The result reads 

r = r [{cos T cos ip — Cq sin T sin ip} p 

+ {C l0 sin T cos ip - (5 4o 5 - C ia Ce cos T) sin ip] q 

+ {S io sin T cos ip + (G 4o 5e + 5 l0 Ce cos T) sin ip} N] 

(7.4) 

where Gi and Si are shorthands for cosio and sinzo, 
respectively. The velocity vector v = dr/dt is given by 



v 



— r [{0 sin sin T sin ip — T(cos cos T sin ip + sin T cos <p) 
—ip(cos sin T cos <p + cos T sin ip) } p 
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+ {6 sin ip(—(d sin 6 cos T + Si cos 6)) + Ci T(cos T cos <p — cos 6 sin T sin ip) 
+y>(cos <p(C io cos 6 cos T — S io sin 6) — C io sin T sin ip)\ q 
+ {6sin</?(Ci cos6 — &,•<, sin6cosT) + Si T (cos T cos ip — cos 9 sin T sin <p) 
+ip(cos >p(Si cos 6 cos T + d g sin 6) — Si sin T sin ip) } AT] . 



(7.5) 



Having inserted above equations into (7.3 1, the final expressions for h x and h + with time dependent O and the case 
of quasi-circular orbits are given by 

h x In -0 ' = i — [(cos T cos ip — cos G sin T sin tp)(Ci ( cos O cos T sin <p> + sin T cos ip) 



-Sig sin 



6 sin (/?)] 

r 2 [— (O sin O sin T sin ip — T(cos O cos T sin ip + sin T cos ip) — c/?(cos sin T cos ip 
cos T sin ip))] (O sin </?(Cio sin(6) cos T + Si cos 6) + Ci T(cos 6 sin T sin ip 
cos T cos ip) + ip(—Ci cos 6 cos T cos(<^) + Si sin 6 cos ip 

C l0 sinTsini/?)) j , 



(7.6) 



/i + k = < (cos T cos ip — cos (3 sin T sin ip) — (sin ?o sm H sin y 

w c 4 it ( r 

— Ci (cos 6 cos T sin ip + sin T cos ifi)) 2 ] 

— r 2 [6 sin ip(Ci sin 6 cos T + Si cos 6) + Ci T(cos 6 sin T sm(tp) — cos T cos ip) 
+ip[—Ci cos cos T cos ip + sin iq sin cos ^ + Ci sin T sin ip)] 2 
+r 2 [0(— sin 0) sin T sin <^ + T(cos cos T sin ip + sin T cos ip) 

+ip(cos sin T cos ip + cos T sin ip)] 2 1 . (7.7) 

I 



VIII. CONCLUSIONS AND OUTLOOK 

In this paper, the EOM of the spins and the orbital 
phase for the conservative 2PN accurate point particle 
were solved for the case of quasi-circular orbits, including 
the leading order spin-orbit interaction. The associated 
gravitational waveforms, h + and h x , in the quadrupo- 
lar restriction are given in analytic form. The spins are 
characterized by their constant magnitudes and 3 essen- 
tial dynamic configuration angles, whose first order time 
derivatives were computed with the aid of Poisson brack- 
ets, and appear to decouple from the orbital phase. Al- 
though these equations are quite complicated and have 
to be integrated numerically in general, they reduce to 
quite simple ones in the case of equal masses and are then 
able to be solved exactly. 

For small deviations from equal masses, a simple per- 
turbative reduction scheme for the EOM can be em- 
ployed. The associated first order corrections to the un- 
perturbed equal-mass solution have been derived. The 
reliability of this solution naturally depends on the pre- 
cision of measurement. The corrections are of the 
same PN order as the unperturbed solutions, multiplied 
by a factor of F = (<5i — S 2 )/(Si + S 2 ). If we set 
W2 = mi (1 + a), we obtain following representa- 



tive pairs (a,F(a)): (0.1, 0.04), (0.2, 0.08), (0.5, 0.17), 
(1.0, 0.29), to give an estimate of the magnitude of the 
perturbation. For the case of a < 0.2, this is below 10 
%, in other cases second-order perturbations may be re- 
quired. 

For a more complete representation, it will be highly 
demanded to include eccentricity (in progress) as well as 
higher-order spin dynamics, which have been found re- 
cently by Steinhoff et al. [15j[T6]. Another important fact 
is that the radiation reaction (RR) has been neglected for 
this analysis. It will be a task to include the energy and 
angular momentum loss due to RR, which is under in- 
vestigation. 
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APPENDIX A: SOLUTION OF THE FULL EOM 
BY LIE SERIES 



For a consistency check, let us solve the problem per- 
turbatively using the Lie series formalism |24j and com- 
pare the results with the computation in the previous 
section. The idea is to associate a linear differential op- 
erator I? to a system of differential equations and to apply 
this operator in an exponential series to the initial val- 
ues. Successive computing of the addends will give the 
perturbed special solution to the required order. Let us 
suppose T> to have the explicit form 



T> = ai(x) 



d 



a n (x) 



d 

dx n 



(Al) 



where n is the number of independent variables and x = 
{xi} with i — (l,...,n). The oti are functions of these 
variables. Then the operator T>, applied to the variable 
Xi, will give 



Vx l = ai(x) . 



(A2) 



Under certain assumptions (holomorphy of the ai), the 
scries 



' ' " - E -jrK*) = f^+tvfix^-v^fix). 

i/=0 



(A3) 

converges absolutely and uniformly for some \t\ < T. 
Defining X to be, in components, 



Xi = (e Xi) with Xi\ t =o = x% , 
the following "exchange relation" 

F(X) = F{e tv x) = e tv F(x) 



(A4) 



(A5) 



holds for the region of convergence. Computing the time 
derivative of the elements Xi, one can use the latter re- 
lation for the operator dt as the function F and obtains 



dXj 
dt 



e tp [Vx i ]=e tp [a i {x)] = a i {X). 



(A6) 



This shows that the are solutions to Eq. . ( A6 ) in the 
region of convergence for the time t. 

The next step is to split the operator D into one part 
T>i, of which the solutions are exactly known, and an- 
other part T>2 perturbing this system of differential equa- 
tions, both supposed to be holomorpic functions in the 
same surrounding of the point x — {xi, ...x n }. Let us 
define the solution to the operators T>\ and T> 2 as 



X<°>(t,x) 
X(t,x) 



e tv *x. 



e tv x 



(A7) 
(A8) 



Inside the region of convergence, the series for X can be 
resummed arbitrarily and cast into another more useful 
form 

X(t,x) =X {0) (t,x) 

v=0 J ° V - 

The series runs over the label v and the operator T>2D V 
has to be applied to x first, the unperturbed solution 
g(t, x) to be inserted afterwards, see [24 for further in- 
formation. The operator V 2 can be shifted before the 
summation, which itself can also be exchanged with the 
integration, and what remains is 

X{t,x) = x(°\t,x)+ I drV 2 X{t-r,x)\ x=xm ^ x) , 
Jo 

(A10) 

This is an integral relation which can be solved itera- 
tively. To any order, for example, the solution reads 

X {1 \t,x) = X {0 \t,x) 

' [D 2 X^(t-T,x)} x=x( ou T . x) }dT, 



I 

Jo 



X {u+1 \t,x) = X {0) (t,x) 

+ I [D 2 X^(t - T,x)] x=x( o HT . x) ]dT . 

Jo 

(All) 

For convergence issues, we note that this expression con- 
verges at least where the double series, Eq. (A8l, con- 



verges absolutely [25]. For a satisfying application of this 
algorithm, the operator T>2 hast to be small; that means 
that the functions a\ (x) (the superscript 2 stands for 
the association to the second operator) are smaller in 

their magnitude in comparison to the coefficients ap{x). 

This algorithm applies excellently to the problem 
of a binary of arbitrarily configurated spins with un- 
equal mass distribution, slightly deviating from the exact 
equal-mass case. The latter is already solved in [TB], and 
what remains is to include perturbations. We will, for the 
time being, resort to the first order of the approximation 
scheme ( |All I to give a representative computation. Of 
course, the results will not sufficiently reflect the physics 
of the system after a long elapsed time and has to be 
expanded for further investigations. 

To first-order in spin-orbit interactions, the motion of 
the spinning binary can be split into the equal mass 
spin-orbit evolution completed by the remainder built 
from the difference in the masses. Let us choose X = 
{O, T, <j> s } as the functions to be evolved, then the EOM 
for X, after the split, symbolically read 



o = xiT 1 (e,0 s ) + x 2 r 2 (o,0 s ) 



(A12a) 
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T = xi«i(e > &) + X2«a(e,0.) > (A12b) 

0. = xi^i(e,^) + X2^2(e,^). (Ai2c) 

The operators 2?i and X>2, therefore, read 

2?! = Xl (Tide + U x d T + Vidt.) , (A13) 

T>2 = X2 (T 2 d e + U 2 dy + V 2 d^) . (A14) 



For the full motion, Eqs. ( A12 1, then X(t) is given by 
the Lie series 

X(t) = e* {Vl+V2) X(t = 0) = e* Pi+^x . (A15) 



The relation for the perturbative functions can be 
computed using the unperturbed one, associated with 
the equal-mass case. The generic angles therein, 
X (0) = {e(°),T(°),0i O) }, read 



T<°>(t) = n T t +T , 

4°\t) = a^t + <f> s0: 
e<°>(t) = e , 



(A16) 
(A17) 
(A18) 



with constant angular velocities, given by Eqs. (4.9 1. The 
first order solutions formally read 



0«(t)-0(°)(t)= / {V 2 e^(t-r,x)} dr, (A19a) 

T«(t) - T<°>(t) = f o {V 2 T«9(t - r, *)}_ x(0)(t ^ ^ , (A19b) 

^ 1 )(i)-^°)(i)= /'{^^(i-r,*)} dr. (A19c) 

All perturbing functions, computed by Eqs. ( A19 1, are in complete agreement with the ones in section [v] 
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